Estimates for the Spectral Condition Number of Cardinal 
B-Spline Collocation Matrices (Long version)* 

Vedran Novakovicj Sanja Singer* and Sasa Singer^ 
April 26, 2010 



Abstract 

The famous de Boor conjecture states that the condition of the polynomial 
B-spline collocation matrix at the knot averages is bounded independently of the 
knot sequence, i.e., depends only on the spline degree. 

For highly nonuniform knot meshes, like geometric meshes, the conjecture is 
known to be false. As an effort towards finding an answer for uniform meshes, 
we investigate the spectral condition number of cardinal B-spline collocation 
matrices. Numerical testing strongly suggests that the conjecture is true for 
cardinal B-splines. 
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1 Introduction 

We consider the classical Lagrange function interpolation problem in the following 
"discrete" setting. Let tx,...,t v be a given set of mutually distinct interpolation 
nodes, and let fx, . . . , /„ be a given set of "basis" functions. 

For any given function g we seek a linear combination of the basis functions that 
interpolates g at all interpolation nodes, 

V 

^Vjfjin) = g(n), i = i,...,u. 

3=1 

The coefficients yj can be computed by solving the linear system 

Ay = g, 
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where A is the so-called collocation matrix containing the values of the basis functions 
at the interpolation nodes 

If the basis functions are linearly independent on {n, . . . ,t u }, the matrix A is non- 
singular, and the interpolation problem has a unique solution for all functions g. 

The sensitivity of the solution is then determined by the condition number k p (A) 
of the collocation matrix 

Kp(A) = ||^4||p||j4 (1) 

where || |L denotes the standard operator p-norm, with 1 < p < oo. 

The polynomial B-splines of a fixed degree d are frequently used as the basis 
functions in practice. Such a basis is uniquely determined by a given multiset of 
knots that defines the local smoothness of the basis functions. The corresponding 
collocation matrix is always totally nonnegative (see [3 [15]), regardless of the choice 
of the interpolation nodes. 

A particular choice of nodes is of special interest, both in theory and in practice, 
for shape preserving approximation. When the nodes are located at the so-called 
Greville sites, i.e., at the knot averages, the interpolant has the variation diminishning 
property. Moreover, for low order B-spline interpolation, it can be shown that ^(A) 
is bounded independently of the knot sequence [6] . 

On these grounds, in 1975, Carl de Boor [S] conjectured that the interpolation 
by B-splines of degree d at knot averages is bounded by a function that depends 
only on d, regardless of the knots themselves. In our terms, the conjecture says that 
Hoo(A) or, equivalently, H^^Hoo is bounded by a function of d only. The conjecture 
was disproved by Rong-Qing Jia [TT] in 1988. He proved that, for geometric meshes, 
the condition number Koo(A) is not bounded independently of the knot sequence, for 
degrees d > 19. 

Therefore, it is a natural question whether exists any class of meshes for which de 
Boor's conjecture is valid. Since geometric meshes are highly nonuniform, the most 
likely candidates for the validity are uniform meshes. 

Here we discuss the problem of interpolation at knot averages by B-splines with 
equidistant simple knots. The corresponding B-splines are symmetric on the support, 
and have the highest possible smoothness. It is easy to see that the condition of this 
interpolation does not depend on the knot spacing h, and we can take h = 1. So, 
just for simplicity, we shall consider only the cardinal B-splines, i.e., B-splines with 
simple knots placed at successive integers. It should be stressed that the only free 
parameters in this problem are the degree d of the B-splines and the size v of the 
interpolation problem. Our aim is to prove that the condition of A can be bounded 
independently of its order. 

The corresponding collocation matrices A are symmetric, positive definite, and 
most importantly, Toplitz. But, it is not easy to compute the elements of A -1 , or 
even reasonably sharp estimates of their magnitudes. So, the natural choice of norm 
in (HD is the spectral norm 

k 2 (A) = ( 2 ) 

Omm\sl) 

where c max (A), and 

(J min( j 4) denote the largest and the smallest singular value of A, 

respectively. 
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For low spline degrees d < 6, the collocation matrices are also strictly diagonally 
dominant, and it is easy to bound K2(A) by a constant. Consequently, de Boor's 
conjecture is valid for d < 6. 

For higher degrees, the condition number can be estimated by embedding the 
Toplitz matrix A into circulant matrices of higher orders. The main advantage of this 
technique, developed by Davis [4] and Arbenz [2], lies in the fact that the eigenvalues 
of a circulant matrix are easily computable. The final bounds for K2(A) are obtained 
by using the Cauchy interlace theorem for singular values (see [TO] for details), to 
bound both singular values in ([2]). 

The paper is organized as follows. In Section 2 we briefly review some basic 
properties of cardinal B-splines. The proof of de Boor's conjecture for low degree (d < 
6) cardinal B-splines is given in Section 3. In Section 4 we describe the embedding 
technique and derive the estimates for K2(A). 

Despite all efforts, we are unable to prove de Boor's conjecture in this, quite prob- 
ably, the easiest case. The final section contains the results of numerical testing that 
strongly support the validity of the conjecture, as well as some additional conjectures 
based on these test results. 

2 Properties of cardinal splines 

Let Xi = xq + ih, for i = 0, . . . , n, be a sequence of simple uniformly spaced knots. 
This sequence determines a unique sequence of normalized B-splines JVg , . . . , N^^^ 
of degree d, such that the spline Nf is non-trivial only on the interval (xj, Xi + d+i)- 

Each of these B-splines can obtained, by translation and scaling, from the basic 
B-spline Q d with knots i = 0, . . . , d + 1, 



Here, (x — i)+ denotes the truncated powers (x — i)+ = (x — i) d (x — for d > 0, 
while 






x < i 



x > i. 



The normalized version of the basic spline is defined as 



N d (x) = (d+l)Q d (x). 
From and (j3|), we obtain the normalized B-spline basis {Nf}: 



(4) 




If the interpolation nodes T{ are located at the knot averages, i.e. 




x i+ i j h x i+d 

d 



Xi + h 



d+1 
2 



i = 



. . . , n — d — 1, 



(5) 



then 



'% < x i < x i+d+l> 
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and the Schonberg-Whitney theorem [8J guarantees that the collocation matrix is 
nonsingular. Moreover, this matrix is totally nonnegative [T2|, i.e., all of its minors 
are nonnegative. Due to symmetry of B-splines on uniform meshes, the collocation 
matrices are also symmetric and Toplitz. So, we can conclude that cardinal B-spline 
collocation matrix is Toplitz, symmetric and positive definite. 

It is easy to show that the elements of the collocation matrix do not depend on 
the step-size h of the uniform mesh, so we take the simplest one with h = 1 and 
Xi = i. Such B-splines are called cardinal. The interpolation nodes ([5]) are integers 
for odd degrees, while for even degrees, the interpolation nodes are in the middle of 
the two neighbouring knots of the cardinal B-spline. 



Nfl\) <V/(x) 




1234567 12345678 



Figure 1: The first four cardinal splines Nf of degree d = 3 and 4, respectively. Big 
black dots denote the spline values at the knot averages. 

The normalized basic cardinal spline N d suffices to determine all basis function 
values at the interpolation nodes 

N d (x*) = N d (x*_ i ). 

The general de Boor-Cox reccurence relation [5J, written in terms of the degree 
of a spline is: 

d _ xN^ 1 (x) + (d + 1 - x)N d - 1 (x - 1) 



N a (x) 



d 



(6) 



Note that the elements of a collocation matrix are rational, because the interpolation 
nodes are rational, and the de Boor-Cox recurrence formula (JBJ) involves only basic 
arithmetic operations on rational coefficients. These elements are therefore exactly 
computable in (arbitrary precise) rational arithmetic. 

3 Low degree cardinal B-splines 

Let tf = N d (x*) denote the values of cardinal B-splines at knot averages (see (JSJ). 
Then, the cardinal B-spline collocation matrix A with interpolation nodes x* is a 
banded Toplitz matrix of order n — d, to be denoted by 

\ ••• ti o 



T 



'o<. 



(7) 



The matrix T d is represented by its first row, usually called the symbol, 



t=(t d ,...,t d r ,0,...,0), t€ 



nn—d 
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It is useful to note that each B-spline of degree d > is a unimodal function, i.e., 
it has only one local maximum on the support. In the case of cardinal B-splines, we 
have already concluded that the splines are symmetric, and therefore the maximum 
values of N d is attained at the middle of the support, for x = (d+l)/2. The maximum 
value is 

N d (^) = N*( X * ) = 4. 

Furthermore, unimodality implies that the values of the spline N d are decreasing in 
the interval [(d + l)/2, d + 1], so 

4 > 4 > • • • > *r- (8) 

To esitmate the condition number of a cardinal B-spline we need to bound both 
the minimal and the maximal singular value of T d . For a symmetric and positive 
definite matrix, the singular values are eigenvalues. Therefore, the bounds for the 
eigenvalues of T d are sought for. From the Gersgorin bound for the eigenvalues, and 
the partition of unity of the B-spline basis, we obtain an upper bound for A max (T^) 

A max (T n d ) <4 + 2(tf + ...+4) = l, (9) 

Similarly, we also obtain a lower bound for X m i n (T d ), 

Amin(T n d ) > 4 - 2(tf + • • • + 4) = 24 - 1, (10) 

which is sensible only if T d is strictly diagonally dominant. Strict diagonal dominance 
is achieved only for B-spline degrees d = 1, ... ,6 (easily verifiable by a computer). 
The corresponding Gersgorin bounds are presented in Table [TJ This directly proves 
de Boor's conjecture for low order B-splines. 



n\d 


2 


3 


4 


5 


6 


64 


1.998136 


2.994873 


4.785918 


7.466648 


11.727897 


128 


1.999541 


2.998757 


4.796641 


7.492176 


11.785901 


256 


1.999886 


2.999694 


4.799180 


7.498105 


11.799106 


512 


1.999971 


2.999924 


4.799797 


7.499534 


11.802256 


1024 


1.999993 


2.999981 


4.799950 


7.499884 


11.803026 


2048 


1.999998 


2.999995 


4.799987 


7.499971 


11.803216 


GB(d) 


2 


3 


f| ps 5.052632 


10 


^ ps 45.354331 



Table 1: Comparison of the actual condition numbers K2(T d ), for d = 2, . . . , 6, n = 
64, ...,2048, and the bounds GB(d) for K2(T d ), obtained by the Gersgorin circle 
theorem. 

Note that in the case of tridiagonal Toplitz matrices, i.e. for d = 2,3, and, thus, 
r = 1 in j7|), the exact eigenvalues are also known (see Bottcher-Grudsky [3]) 

Afc(T^) = ^ + 2ifcos ^— , d = 2,3, k = l,...,n-d. 
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The largest and the smallest eigenvalue can then be uniformly bounded by 

Amax(T n d ) = t d + 2tf COS < tg + 2tj, 

Amin(T„ d ) = 4- 2t( COS n _ 7T d+1 >4~ 2 4 > 0, 

These uniform bounds are somewhat better than those obtained by the Gersgorin 
circles. 

4 Embeddings of Toplitz matrices into circulants 

When the degree of a cardinal B-spline is at least 7, the eigenvalue bounds for Toplitz 
matrices can be computed by circulant embeddings. First, we will introduce the 
smallest possible circulant embedding, and give its properties. Then we will present 
some other known embeddings, with positive semidefinite circulants. 

To obtain a bound for A m i n (T^), the collocation matrix T d is to be embedded into 
a circulant 
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m = n — d + r. (11) 



It is obviously a Toplitz matrix with the following symbol 

t = (t , . . . ,t r , 0, . . . , 0, t r , . . . ,ti), t£R n d+r . 

This circulant is called a periodization of T d by Bottcher and Grudsky [3]. 

The bounds (j9|)- (fT0|) for the eigenvalues of T d are also valid for C^. Moreover, 
is doubly stochastic, always having A max (C^J = 1 as its largest eigenvalue. In- 
terestingly enough, the upper bound ([9|) is attained here (the Gersgorin bounds are 
rarely so sharp). 

The symmetry of T d immediately implies the symmetry of C^, and we can con- 
clude that the eigenvalues of are real, but not necessarily positive. For symmetric 
matrices, the singular values are, up to a sign, equal to the eigenvalues, so 

a i {C d m ) = \\ l {C d m )\. (12) 

If the eigenvalues of the circulant C^ are known, the spectrum of embedded T d can 
be bounded by the Cauchy interlace theorem for singular values, applied to C^. 

Theorem 4.1 (Cauchy interlace theorem). Let C G c mxn fr e given, and let Cg denote 
a submatrix of C obtained by deleting a total of £ rows and/or £ columns of C . Then 

o-fc(C) > a k (Ci) > a k+ i(C), k = 1, . . . ,min{m,n}, 

where we set o-j(C) = if j > min{m,n}. 
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The proof can be found, for example, in [10, page 149]. 

If we delete the last r rows and columns of C^, we obtain Tf t . The Cauchy 
interlace theorem will then give useful bounds for a m - m (T^) = A m i n (T^), provided 
that C^li is nonsingular. Moreover, if we delete more than r last rows and columns of 
C^j, we obtain bounds for Toplitz matrices Tj}, of order k — d, for k < n, 



„ Crd\ a ma,x{Tj}) 0"max(C^j) 
«2UfcJ = : < 



&rain(Tj}) 0"min(C m ) minj \Xj(C m 



(13) 



Now we need to calculate the smallest singular value of and show that it is 
non-zero. 

The eigendecomposition of a circulant matrix is well-known (see [U[2])- A circu- 
lant C of order m, defined by the symbol (cq, . . . , c m _i), can be written as 



m— 1 



where 



3=0 
1 



n 



i 



The spectral decomposition of II is IT = FQF* , where 

2tti 



O = diag(l, u), u) ,...,ui 



m— 1 \ 



m 



while 



F, 



3,k 



m 



uj kj , 0<k,j <m-l. 



Hence, C can be decomposed as 



m— 1 



C = FAF*, A = diag(A , . . . , \ m -i) = ^ CjW . 
The eigenvalues of a real symmetric circulant C are real, and given by 

m— 1 



Afc(C) = co + c j cos 

3=1 



2nkj 



, k = 0, . . . , m — 1. 



(14) 



They can also be viewed as the discrete Fourier transform (DFT) of the symbol 
[co, • • • , c m _i). 

For real and symmetric C, i.e., when ct = c m -k, for k = 1, . . . , m — 1, from ([Tj 
it also follows that 



m— 1 ^ j^, • m— 1 

Afc(C) = C + ^ Cj COS = Cq + ^ 

i=i ; " j=i 



2irk(m — j) 



m 



Am-fe(C)- 



So, all the eigenvalues, except Aq(C), and possibly Am(C), for even to, are multiple. 
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Therefore, the eigenvalues of the circulant from (|lip are 

X k (Ci) = 4 + 2 Vtfcos^^, fc = 0,...,m-l. (15) 
i=i 

For prime orders m, the nonsingularity of is a consequence of the following 
theorem from [§]. 

Theorem 4.2 (Geller, Kra, Popescu and Simanca). Lei m be a prime number. As- 
sume that the circulant C of order m has entries in Q. Then det C = if and only 
*f 

m—l 

X = ^2 a = 0, 

3=0 

or all the symbol entries q are equal. 

If m is prime, then we must have detC^ ^ 0, since ([SD implies that Cj's are not 
equal, and from (fT5|) we get 

r 

Ao = c + 2^ Ci = l/0. 
j=i 

Theorem 14.21 suggests how to get the nonsingular embedding of T!~. First, should 
be embedded into the Toplitz matrix Tp, of order p — d, where p > n is chosen so 
that m = p — d + r is a prime number. Then, is embedded into the circulant C^. 

The other possibility is to embed into the smallest circulant matrix C^, as in 
(fTT|) . and calculate its eigenvalues from (fT5|) . in hope that is nonsingular. In this 
case, extensive numerical testing suggests that is always positive definite, but we 
have not been able to prove it. 

There are also several other possible embeddings that guarantee the positive 
semidefiniteness of the circulant matrix C. 

The first one, constructed by Dembo, Mallows and Shepp in [7|, ensures that 
the positive definite Toplitz matrix T, of order n, can be embedded in the positive 
semidefinite circulant C, of order m, where 

m>2^n + K 2 (T)^= S J. (16) 

A few years later, Newsam and Dietrich [14j reduced the size of the embedding to 



m > 2W6n2 + ^{rf—^ . (17) 

Note that among all positive semidefinite matrices C of order greater or equal m, we 
can choose one of prime order. This embedding will be positive definite according to 
Theorem 14.21 It is obvious that embeddings (fT6|) ~ ([T7|l are bounded by a function of 
the condition number of T, i.e., the quantity which we are trying to bound. 

Ferreira in [S] embeds a Toplitz matrix T of order n, defined by the symbol 
t = (to, . . . , t r , 0, . . . , 0) G W 1 , into the circulant C of order m = 2n, 



C 



T S 
S T 



(18) 
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where the symbol of the Toplitz matrix S is s = (0, . . . , 0, t r , . . . , t\) G M. n . 

If we take T = from ([7D, the only difference between embeddings and (|18|) 
is in exactly n — d — r zero diagonals, added as the first diagonals of S. A sufficient 
condition for positive semidefiniteness of C is given by the next result. 

Theorem 4.3 (Ferreira). Let C be defined as in A18\) , and let b T = [to, . . . , t n _ i], 
c T = [t n -i, ■ ■ ■ ,t\]. If T is positive definite, and \b T T~ 1 c\ < 1, then C is positive 
semidefinite. 

Once again, there is no obvious efficient way to verify whether the condition 
\b T T~ 1 c\ < 1 is fulfilled or not. 



5 Conjecture about the minimal eigenvalues 

Extensive numerical testing has been conducted, by using Mathematica 7 from Wol- 
fram Research, for the symbolic, arbitrary-precision rational, and machine-precision 
floating-point computations. Whenever feasible, the full accuracy was maintained. 
Owing mostly to the elegance and the accuracy of these results, insight into and the 
following conjecture about the spectral properties of the collocation matrices and the 
corresponding periodizations were obtained. 

Conjecture 5.1 (The smallest eigenvalue of a circulant). The circulant from 
< T77]) is always positive definite, and the index \i of its smallest eigenvalue A M (C^J is 
always the integer nearest to m/2, i.e., 



•V(C m ) 



i=i 



A m ±i (C%) =4 + 2^(-iy'^cos ( , m odd, 



A f (0=^ + 2^(-l)^, 

3=1 



(19) 



m even. 



Figure [2] illustrates both cases of Conjecture 15.11 



3 4 5 6 7 8 9 10 



13 14 15 16 17 



0.4 - 
0.2 - 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 



Figure 2: The eigenvalues (black dots) Afc(C^) for spline of degree d = 7 with n = 
23, 24, respectively. The associated circulants have order m = n — 7 + 3, i.e., 19 and 
20. Note that for m = 20 there is only one minimal eigenvalue, while for m = 19 we 
have two minimal eigenvalues. 

For even m, A /Lt (C~ l ) (and, therefore, /^(C-J) depends solely on d, i.e., the order 
m of a circulant is irrelevant here. Moreover, for m odd and even alike, the limiting 
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value of A^(C~J is the same: 



lim X tl (Ci) = t d + 2j2(-m 



(20) 



Hence, the notation AJ^ is justified, since that value is determined uniquely by the 
degree d of the chosen cardinal splines. This is consistent with de Boor's conjecture. 

The equations (|19j) and (|20[) provide us with efficiently and exactly computable 
estimates of the spectral condition numbers of large collocation matrices T d . As 
demonstrated in Figure [3] and Table [21 the smallest eigenvalues of the collocation 
matrices converge rapidly and monotonically to the smallest eigenvalues of the cor- 
responding circulant periodizations C^, as well as to the limiting value ([20]) . 




320 512 704 896 1088 1280 1472 1664 1856 2048 



Figure 3: Spectral condition numbers of Toplitz matrices (lower, brighter line), 
and the circulant periodizations (solid black line). The constant function denotes 

I/A" 



It is worth noting that the spectral bounds obtained in such a way for lower 
degrees (d = 2, . . . , 6) of cardinal B-splines are quite sharper than those established 
by the Gersgorin circle theorem (cf. Tableland Table [2]), at no additional cost. 

Since tfj are rational numbers, (|20p is useful for the exact computation of AJ^,. But, 
in floating-point arithmetic, the direct computation of A^, from (|20p is numerically 
unstable, as it certainly leads to severe cancellation. 

It can be easily shown from (|3|) or (|6|) that the smallest non-zero value of the 
cardinal B-spline of degree d at an interpolation node is: 



'N d (l) = for odd d, 

iVd (l) = 2^I' for even d. 



Moreover, all other values t| in ([20]) and, consequently, A^ are integer multiples of t d . 
With that in mind, yet another, somewhat surprising conjecture emerged from the 
test results: 

(t d -T d = ±T d , dodd, 

A - = , * i (21) 
\t d -2 d E d = ±E d , deven, 



10 



n\d 


T 2 


r 2 


± n 


r 5 


rp6 

± n 




64 


1.998137 


1.998758 


7.466648 


7.472749 


11.72790 


11.74214 


128 


1.999541 


1.999694 


7.492176 


7.493492 


11.78590 


11.78866 


256 


1.999886 


1.999924 


7.498105 


7.498410 


11.79911 


11.79971 


512 


1.999971 


1.999981 


7.499534 


7.499607 


11.80226 


11.80240 


1024 


1.999993 


1.999995 


7.499884 


7.499902 


11.80303 


11.80306 


2048 


1.999998 


1.999999 


7.499971 


7.499975 


11.80322 


11.80322 






2.000000 




7.500000 




11.80328 


n\d 


n 




7^21 

n 




7^30 

n 


/^30 


64 


45.04067 


45.17179 


9012.21 


9543.49 


371000.6 


502472.1 


128 


45.57648 


45.59721 


10100.96 


10150.47 


569223.5 


579852.3 


256 


45.69092 


45.69486 


10273.67 


10279.58 


594976.6 


596037.0 


512 


45.71737 


45.71822 


10308.14 


10309.00 


599497.1 


599628.0 


1024 


45.72373 


45.72393 


10315.86 


10316.01 


600450.4 


600469.5 


2048 


45.72529 


45.72534 


10317.69 


10317.72 


600669.7 


600673.0 






45.72581 




10318.28 




600739.5 



Table 2: Comparison of the spectral condition numbers ^(T^) and ^(C^J, for 
d = 2, 5, 6, 9, 21, 30, n = 64, . . . , 2048, m = n - d + r, and 1/A^. 



where, as in |13| . T n are the tangent numbers, and E n are the Euler numbers, defined 
by the Taylor expansions of tani and seci, respectively, 



taut = V]r n — , 

n=0 



see 



nl 

n=0 



These numbers are also related to the sequences A000182 (the tangent or "zag" 
numbers), A000364 (the Euler or "zig" numbers) and A002436, from [16J. 

If true, ()2ip would be of significant practical merit, for there exist very stable and 
elegant algorithms for calculation of T n and E n by Knuth and Buckholtz |13j . So, it 
deserved an effort to find the proof. 

A unifying framework for handling both cases is provided by the Euler polyno- 
mials E n (x), defined by the following exponential generating function (see [U 23.1.1, 
p. 804]) 



o xt °° 



e + 



n=0 



nl 



(22) 



which is valid for \t\ < it. 

First, note that T 2 k = E 2 k+i = 0, for all k > 0. The remaining nontrivial values 
can be expressed in terms of special values of Euler polynomials. For the tangent 
numbers, we have 

T 2k+1 = (-l) k 2 2k+1 E 2k+1 (l), k>0. (23) 
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This follows easily, by comparing the Taylor expansion of 1 + tanh t 



1 + tanh t 



2e 



e 2t + 



2t 00 



•'71 



t 



2k+l 



2k+l 



k=0 



(2k + l)\ 



and ()22p . with x = 1 and 2£, instead of t. Similarly, by comparing the Taylor 
expansion of sech t 



secht 



2e l 



e 2t + 



- = ^(-l)^ 2fc 



fc=0 



i-2 k 



(2k)\ 



and ([52]). with x = 1/2 and 2i, instead of i, we get 



E 2k = (-l) h 2 2k E 2k 



k > 0. 



The following identities will also be needed in the proof of ([21 
Lemma 5.2. Lei d> be a non-negative integer. Then 



= 0, 



d+1 



E(- 1 )'( d t 1 )^+ 1 ) = ' 



1=0 



for all n = 0, . . . , d. 

Proof. Consider the function g^ defined by 



9d(t) :-- 



2(1 - e 



t\d+l d +! 



e* + 1 



£=0 



d + 2e 



.ft 



e* + 1 



From (|22p with x = £, the Taylor expansion of can be written as 



9d 



w = E 



n=0 



f><r) 



4W 



so 



£=0 
d+1 



n! 



B-wWU-EC-D'fT) 

e=o \ 1 / 



E n (£), n > 0. 



On the other hand, the Leibniz rule gives 



m=0 



e* + 1 



If n < cZ, then D m [(1 — e*) d+1 ] is always divisible by (1 — e l ). Hence, 



D n g d {t) 

which proves the first identity ([251 



t=o 



0, n = 0, . . . , d, 



(24) 



(25) 
(26) 
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The second one follows similarly, by considering 

h d (t) := g d (t) - g d+1 (t) = K > = J^(-l) M J 

The Taylor expansion of h d is then given by 



e* + 1 



M*) = E 



n=0 



i=0 



??■! 



If n < d, from the first part of the proof, it follows immediately that 

D n h d (t)\ t=0 = D n g d (t)\ t=Q - D n g d+1 (t)\ t=Q = 0, 
which proves (|26|) . □ 

Finally, we are ready to prove the conjecture (|21|) . 

Theorem 5.3 (Relation to integer sequences). The following holds for all cardinal 
B-spline degrees d > 

^ d _ 1_ \T d , d odd, 
°° d\ yE d , d even. 

Proof. To simplify the notation, let L d := d\ A^,. Due to the symmetry of interpola- 
tion nodes, the sum in (1201) can be written as 



j = _ r . . . r , 

where r = |_rf/2j . From ([3]) and (jJJ), it follows that 

d + l /j , i\ / j , 1 \ rf 



*-*t<-<r)(>+*H: 



Then 

j=-r £=0 

Let d be odd, d = 2k + 1, with fc > 0. Then r = k and (d+ l)/2 = k + 1, so ([2 
becomes 

j=-k e=o ^ ' 

From the definition of truncated powers with positive exponents, the second sum 
contains only the terms with j — £ + k + 1 > 0, i.e., for I < j + k. By changing the 
order of summation, we get 

2k k 

w=£(-i/r+ 2 ) E (-iyu-i+k+i) 2k+1 . 

e=o ^ ' j=£-k 
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Then we shift j by k — £ + 1 , so that j starts at 1 , to obtain 

2k /ri , , _ x 2k+l- 



(2fe+l-£)-j -2fe+l 



£=0 v ' j=l 

The second sum can be simplified as (see [U 23.1.4, p. 804]) 

2k-\-l— t 

E (-1)^+1-^)^/^1 = I (^ 2fc+1 (2fc + 2 - £) + (-l) 2fe+2 -^ 2fc+1 (i; 
Hence 



i 



( 2 \ + 2 ) we* + 2 - o + iwu) E ( 2 "+ 2 ) 



.£=0 v ' 1=0 

By reversing the summation, from (|25|) with d = 2A; + 1 and n = d, we conclude that 

2k /r , 7 , „ N 2fc+2 



£(-tf ( 2 \ + 2 ) we* + 2 - = E f 2 \ + 2 ) WW 

£=0 V J 1=2 V 7 

B-i)'( 2 \ +2 WiW- 

0— n \ / 



i=o 

Since i? 2 fc+i(0) = — E 2 fe+i(l), by using ((23j), we have 

^2fc+l = - ^ - E ' 2fc + 1 ( 1 ) E ( / ) = ( _1 ) fc22fc+l£ ' 2fc + 1 ( 1 ) =7 2fc+l- 

This proves the claim for odd values of d. 

Let d be even, d = 2k, with > 0. For cZ = 0, it is obvious that Lq = t$ = 1 = Eq, 
so we may assume that k > 0. Then r = k and (cZ + l)/2 = + 1/2, so (|27p becomes 

fc 2fc+l /o; , i\ / i \ 2fc 

j=-k 1=0 v 7 v J + 

The second sum contains only the terms with j — £ + A; + 1/2 > 0, i.e., for I < j + k. 
By exactly the same transformation as before, we arrive at 

2k , , x 2fc+l-£ , . 2fe 

^» = (-i)*E(- 1 ) / ( / ) E (-i) (2fc+1 -^' (j - y • 

Now we expand the last factor in terms of powers of j. Then L 2 & can be written as 
^k = (-l) k f:( 2 ")(-l) 2k ~ n S 2k , n , (28) 

with 



n=0 v 



2k , , v 2fc+l-£ 

s 2M = E(-i)v ^ j ) E (-i) (2fc+1 - £) - j r, n = o,...,2*. 

£=0 ^ ' 7=1 



14 



Like before, the second sum can be simplified as 

2fc+l-£ 

£ = | (JSW(2* + 2 - £) + (-l) 2fc+2 -^ n (l)) , 

which gives 

2fc /„, ,\ 2fc 



5*2A;,n — ^ 



EM)' f% + ') ^ + 2 - £) + E n (l) £ ( 2fc + X ) 

tf— n \ / 0— n V / 



J=0 v ' £=0 

By reversing the summation, from (j26|) with d = 2/c, for n = 0, . . . , d, we see that 



1=0 

Therefore 



1=0 

From (1281) we obtain 



5» >B = i^(l)^ 1 ( 2fc + 1 ) =2*^(1). 



2k /o; \ / i \ 2fc-n 



^=(-d^e( 2 „>.w(4) 

n=0 v 7 v 7 



Finally, by using [Tj 23.1.7, p. 804]) 



Together with (|24p . this gives 

L 2fc = (-l)^ 2fc (i)=^. 

This completes the proof for even values of d. □ 

We would like to conclude with an observation that, to the best of our knowledge, 
scarcely any result could be found about sufficient conditions for the non-negativeness 
of the DFT in terms of its coefficients, apart from the classical result of Young and 
Kolmogorov (cited in Zygmund [TTJ, page 109]): 

Theorem 5.4. For a convex sequence (a n ,n 6 N), where lim a n = 0, the sum 



oo 

.. cos (nx) 

n=l 



converges (save for x = 0), and is non-negative. 

Here, a sequence is convex if A 2 a n > for all n, with Aa n = a n — a n+ \. 

Convexity is not fulfilled in the case of cardinal B-spline coefficients, since there 
is always one inflection point on each slope of the spline. And yet, our numerical 
experiments strongly suggest that the class of series with a positive DFT is worth 
investigating further, for the theoretical and practical reasons alike. 
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